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■ We consider a general, classical theory of gravity in n dimensions, arising from 

I a diffeomorphism invariant Lagrangian. In any such theory, to each vector field. 
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on spacetime one can associate a local symmetry and, hence, a Noether current 



Q I (n — l)-form, j, and (for solutions to the field equations) a Noether charge (n — 2)- 



form, Q, both of which are locally constructed from ^'^ and the the fields appearing 
in the Lagrangian. Assuming only that the theory admits stationary black hole 
solutions with a bifurcate Killing horizon (with bifurcation surface S), and that 
the canonical mass and angular momentum of solutions are well defined at infinity, 
we show that the first law of black hole mechanics always holds for perturbations to 
nearby stationary black hole solutions. The quantity playing the role of black hole 
entropy in this formula is simply 27r times the integral over S of the Noether charge 
(n — 2)-form associated with the horizon Killing field (i.e., the Killing field which 
vanishes on S), normalized so as to have unit surface gravity. Furthermore, we 
show that this black hole entropy always is given by a local geometrical expression 
on the horizon of the black hole. We thereby obtain a natural candidate for the 
entropy of a dynamical black hole in a general theory of gravity. Our results show 
that the validity of the "second law" of black hole mechanics in dynamical evolution 
from an initially stationary black hole to a final stationary state is equivalent to 
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the positivity of a total Noether flux, and thus may be intimately related to the 
positive energy properties of the theory. The relationship between the derivation 
of our formula for black hole entropy and the derivation via "Euclidean methods" 
also is explained. 

PACS #: 04.20.-q, 97.60.Lf 
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One of the most remarkable developments in the theory of black holes in classical 
general relativity was the discovery of a close mathematical analogy between certain 
laws of "black hole mechanics" and the ordinary laws of thermodynamics. When the 
effects of quantum particle creation by black holes were taken into account, this 
analogy was seen to be of a physical nature, and it has given rise to some deep insights 
into phenomena which may be expected to occur in a quantum theory of gravity. 

The original derivation of the laws of black hole mechanics in classical general relativ- 
ity PI used many detailed properties of the Einstein field equations, and, thus, appeared 
to be very special to general relativity. However, recently it has become clear that at 
least some of the laws of classical black hole mechanics hold in a much more general 
context. In particular, it has been shown that a version of the first law of black hole 
mechanics holds in any theory of gravity derivable from a Hamiltonian . (For the cases 
of (1 + l)-dimensional theories of gravity ^ and Lovelock gravity 0, the explicit forms 
of this law have been given.) Furthermore, analogs of all of the classical laws of black 
hole mechanics have been shown to hold in (1 + l)-dimensional theories 0]. 

However, despite the very general nature of the Hamiltonian derivation |^ of the first 
law of black hole mechanics, there remains one unsatisfactory aspect of the status of 
the first law in a general theory of gravity [^: Although the derivation shows that for 
a perturbation of a stationary black hole, a surface integral at the black hole horizon 
(involving the unperturbed metric and its variation) is equal to terms involving the 
variation of mass and angular momentum (and possibly other asymptotic quantities) 
at infinity, the derivation does not show that this surface term at the horizon can be 
expressed as /t/27r (where k denotes the unperturbed surface gravity) times the variation 
of a surface integral of the form S = Jy,F, where F is locally constructed out of the 
metric and other dynamical fields appearing in the theory. It is necessary that the 
horizon surface term be expressible in this form in order to be able to identify a local, 
geometrical quantity, S, as playing the role of the entropy of the black hole. 

The main purpose of this paper is to remedy this deficiency by showing that in 
a general theory of gravity derivable from a Lagrangian, the form of the first law of 
black hole mechanics for perturbations to nearby stationary black holes is such that the 
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surface term at the horizon always takes the form ^SS, where S" is a local geometrical 
quantity, and is equal to 27r times the Noether charge at the horizon of the horizon Killing 
field (normalized so as to have unit surface gravity). The local, geometrical character 
of S suggests a possible generalization of the definition of entropy to dynamical black 
holes. The relationship between black hole entropy and Noether charge also suggests the 
possibility of a general relationship between the validity of the second law of black hole 
mechanics (i.e., increase of black hole entropy) and positive energy properties of a theory. 
An additional byproduct of our analysis is that it will enable us to make contact with 
the "Euclidean derivation" of formulas for black hole entropy, thereby demonstrating 
equivalence of that approach with other approaches - a fact that is not at all easy to see 
by a direct comparison of, say, references @] and 0. Our considerations in this paper 
will be limited to a general analysis of all the above issues; applications to particular 
theories will be given elsewhere 

Before presenting our new derivation of the first law, we comment upon the status of 
other "preliminary laws" of black hole mechanics in a general theory of gravity. We con- 
sider theories defined on an n-dimensional manifold M with dynamical fields consisting 
of a (Lorentzian) spacetime metric, gab, and possibly other matter fields, such that the 
equations of motion of the metric and other fields are derivable from a diffeomorphism in- 
variant Lagrangian. (Our precise assumptions concerning the Lagrangian will be spelled 
out in more detail below.) We assume that a suitable notion of "asymptotic flatness" is 
defined in the theory. The black hole region of an asymptotically fiat spacetime then is 
defined to be the complement of the past of the asymptotic region. In order to begin 
consideration of the classical laws of black hole mechanics, it is necessary that the event 
horizon of a stationary black hole be a Killing horizon, i.e., a null surface to which a 
Killing vector field is normal. This property is known to be true in general relativity 
by a nontrivial argument using the null initial value formulation 0, so it is not obvious 
that it would hold in more general theories of gravity. Nevertheless, this property auto- 
matically holds for all static black holes (since the static Killing field must be normal to 
the event horizon of the black hole), and, hence, it automatically holds for spherically 
symmetric black holes (in the 0{n — l) sense) and, thus, in particular, for all black holes 
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in (1 + 1) -dimensional theories of gravity. 

The surface gravity, k, at any point, p, of a Kilhng horizon Ti. is defined by (see, e.g., 

CVae = (1) 

where ^'^ is the Kilhng field normal to Ti. The zeroth law of black hole mechanics asserts 
that K is constant over the event horizon of a stationary black hole. The proof of this law 
does make direct use of the specific form of the Einstein field equations and, thus. 



does not appear likely to generalize to other theories of gravity [|TT|. Nevertheless, the 
zeroth law trivially holds for spherically symmetric black holes, and, in particular, in all 
(1 + l)-dimensional theories. 

It is worth noting that the validity of the zeroth law is, in essence, equivalent to the 
statement that - apart from the "degenerate case" of vanishing surface gravity - the event 
horizon of a stationary black hole must be of bifurcate type. Namely, it is easily proven 
that a bifurcate Killing horizon must have constant (and nonvanishing) surface gravity. 



whereas it can be shown |T2| that any Killing horizon with constant, nonvanishing surface 



gravity can be locally extended (if necessary) to a bifurcate horizon. 

It also should be noted that in an arbitrary theory of gravity, a black hole with con- 
stant surface gravity will "Hawking radiate" at temperature k,/2tt when quantum particle 
creation effects are taken into account, i.e., the Einstein field equations play no role in 
the derivation of the Hawking effect. Similarly, the theorems of [0 on the uniqueness 
and thermal properties of quantum states on black holes with bifurcate horizons hold in 
an arbitrary theory of gravity. Thus, k/27i always represents the physical temperature 
of a black hole. 

We turn, now, to the presentation of a new derivation of the first law of black hole me- 
chanics for stationary black holes with bifurcate horizons in a general theory of gravity in 
n-dimensions derived from a diffeomorphism invariant Lagrangian. We shall follow closely 
the framework of Lagrangian field theories developed in WM, with one small change: We 



shall view the Lagrangian as an n-form, L, rather than as a scalar density; similarly, other 



tensor densities of will appear here in their dualized version as differential forms. In 



order to define L, it is necessary to introduce a fixed (i.e., "nondynamical" ) derivative 



operator, Va, on spacetime. It also may be necessary to introduce other "non-dynamical, 
background fields", 7, such as the curvature of Va (if Va is non-fiat); we shall assume, 
however, that any such additional fields, 7, are uniquely determined by Va, and that 
7 changes by a diffeomorphism under the change induced in Va by the action of that 
diffeomorphism. At each point p of spacetime, L then is required to be a function of 
the spacetime metric, Qab, (or, alternatively, of a tetrad or soldering form) and finitely 
many of its (symmetrized) derivatives at p, as well as of other matter fields present in the 
theory and their (symmetrized) derivatives at p, and of 7 at p. Note that no restriction 
is placed upon the number of derivatives of the metric or other fields upon which L can 
depend (other than that this number be finite), so "higher derivative" gravity theories 
are included in this framework. 

In order to reduce the number of symbols and indices appearing in formulas, I shall 
use the symbol "0" to denote all of the dynamical fields, including the spacetime metric. 
We shall restrict attention to diffeomorphism invariant theories, by which we mean that 
for any diffeomorphism, %p : M ^ M, we have, 

Mrm = rm (2) 

Note that on the left side of this equation, tp* is not applied to Va or any other non- 
dynamical fields 7 which may appear in L. Equation (0) can be interpreted as stating 
that - although it may be necessary to introduce Va and/or 7 to define L - L actually 
depends only upon the dynamical fields (p. 

Under a first order variation of the dynamical fields, the variation of L can be put in 
the form (see, e.g., |p!4|), 

6L = E6(f) + d@ (3) 

where summation over the dynamical fields (and contraction of their tensor indices with 
corresponding dual tensor indices of E) is understood in the first term on the right 
side of this equation. The {n — l)-form, 0, is locally constructed from (p and 6(j), but 
is determined by eq.(§) only up to addition of a closed (and, hence, exact |jl^) form 
locally constructed from the fields appearing in L; we shall adopt eq. (2.12) of as 
our definition of 0. The symplectic current {n — l)-form, f2, is defined in terms of the 
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variation of by, 

^^(0,5l0,(520) = (5i[e(0,(520)] -52[e(0,5i0)] (4) 

It should be noted that and O will depend upon the choice of Va in sufficiently high 



derivative theories - although they change only by an exact form, i.e., a "surface 
term", under a change of derivative operator - and they need not be diffeomorphism 
invariant in the sense of eq.(0). Furthermore, and f2 will change when an exact form 
is added to L - with the change in fl being given by an exact form - even though such 
a modification of L has no effect upon the equations of motion, E = 0. 

Now, let C,"" be any vector field on M and consider the field variation 50 = C^(j). The 
diffeomorphism invariance of L implies that under this variation, we have, 

6L = C^L = ■ L) (5) 

where here and below, we make frequent use of the general identity 

C^A = ^-dA + ■ A) (6) 

holding for any differential form A and vector field where "•" denotes the contraction 
of a vector field with the first index of a differential form. Equation shows that the 
vector fields on M constitute a collection of infinitesimal local symmetries in the sense 
of Hence, to each we may associate a Noether current (n — l)-form, j, defined by 



j = 0(0,£50)-^L (7) 

so that j is locally constructed out of the fields appearing in L and A standard 
calculation shows that 

d] = -E£50 (8) 

so that j is closed whenever the equations of motion are satisfied. Since j is closed for all 
it follows [|T3| that there exists an {n — 2)-form, Q - locally constructed out of the 



fields appearing in L and - such that when evaluated on solutions to the equations of 
motion, we have, 

j = dQ (9) 



Since j depends linearly on we adopt the explicit algorithm provided by lemma 1 of 



JTSf to uniquely define Q, from which it follows that Q depends on no more than {k — 1) 
derivatives of ^""j where k denotes the highest derivative of any dynamical field occurring 
in L. (Note, however, that Q is unique up to addition of a closed - and, hence, exact 



Ill5| - (n — 2)-form locally constructed from the fields appearing in L and from so 
the integral of Q over any closed {n — 2)-dimensional surface, S, is uniquely defined by 
eq.(^ alone.) We shall refer to Q as the Noether charge {n — 2) -form |T6[ relative to 
and its integral over a closed surface, E, will be referred to as the Noether charge of S 
relative to 

The key identity upon which our derivation of the first law of black hole mechanics 
will be based is obtained by considering the variation of eq.(^ resulting from an arbitrary 
variation, 50, of the dynamical fields off of an arbitrary solution 0. We have, 

<5j = (5[0(0,/:50)]-e-5L (10) 

(Note that is held fixed in this variation, i.e., we require that 5^°" = 0.) However, by 
eq.(|]), we have, 

e • 5L = e ■ + d&] 

= C^&-di^-@) (11) 

where the equations of motion, E = 0, for and the identity (§) were used in the second 
line. Thus, we obtain, 

5} = 5[0(0, £^0)] - £^[0(0, 50)] + rf(e • 0) (12) 

Note that in eq.(|12|), no restrictions have been placed upon 6(f) or 

Our next step is to identify certain "surface terms" appearing in eq.(|l^. First, we 
require Va to be invariant under the diffeomorphisms generated by This requirement 
holds in the usual case where is taken to be a coordinate vector field and Va is taken 
to be the coordinate derivative operator of that coordinate system; it also will hold in 
our main application below where Va will be taken to be the derivative operator of the 
unperturbed metric and is a Killing field of that metric.) In that case, the first two 
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terms on the right side of eq.(|T^) combine to yield 

(5[e(0, £^0)] - £^0(0, 54>)] = O(0, 50, £^0) (13) 
and eq.(|T^) becomes simply, 

5j = f^(0,50,£g0) + rf(e-0) (14) 

When integrated over a Cauchy surface, C of the unperturbed solution, eg . ([T^ corre- 
sponds to eq.(3.22) of [^, but eq.(|l4|) contains vital additional information concerning 



the "surface term", ■ 0), which did not appear in [Ti], since attention there was 
restricted to the case of compact C. Comparison of eg . {\L^) with Hamilton's equations of 
motion shows that if a Hamiltonian, H, corresponding to evolution by exists on phase 
space, then H must satisfy, 

6H = 6j^}- j^d{i-@) (15) 

where, in this eguation, projection of the right side to phase space (in the manner dis- 
cussed in [|14|) should be understood. This shows that apart from the "surface term" 



■ 0), the Noether current, j, acts as a Hamiltonian density. 

We now further restrict attention to the case where 50 satisfies the linearized equations 
of motion, so that both and its variation are solutions. Then we may replace j and 
its variation by rfQ in eqs.(|l^) and (plSj) . It then can be seen immediately from eq.(|T5|) 
that the Hamiltonian - if it exists - is purely a "surface term". In an asymptotically 
fiat spacetime, it is natural to associate the value of the surface contribution to the 
Hamiltonian from infinity with the corresponding "conserved quantity" associated with 

in the manner of |T^. In other words, if the theory admits a suitable definition of 



the "canonical energy", 8, associated with an asymptotic time translation, t", and of 
the "canonical angular momentum" J , associated with an asymptotic rotation, y)", the 
variations of these quantities should be given by the formulas 

58= f (5Q[t] - t • 0) (16) 

J oo 

6J= - [ 5QM (17) 



9 



where the integrals are taken over an [n — 2)-dimensional sphere at infinity and the term 
(f ■ © does not appear in eq.(|l7D because y?" is assumed to be tangent to this sphere. 
Thus, if one can find an (n — l)-form, B, such that 

S [ t B = [ t & (18) 

the canonical energy and angular momentum can be defined by, 

£ = [ (QM-t-B) (19) 

J = - I (20) 

J oo 

Note that £ corresponds to the "ADM mass" of general relativity plus possible additional 
contributions from any long-range matter fields that may be present; see for explicit 
discussion of the case of the Yang-Mills field. Note also that for the Hilbert Lagrangian 
of general relativity, the expressions Q[t] and — Qiv] correspond - up to numerical 
factors - to the Komar expressions for mass and angular momentum. The presence of 
the "extra term" t • B in eq.(0) accounts for why different relative numerical factors 
must be chosen in the Komar formulas for these quantities.) It is, of course, a nontrivial 
condition on a theory that it admit a notion of asymptotic fiatness such that £ and J 
are well defined. In the following, I shall assume that this is the case, and derive the first 
law of black hole mechanics for such a theory. 

Consider, now, a stationary black hole solution with a bifurcate Killing horizon, with 
bifurcation {n — 2)-surface S. Choose to be the Killing field which vanishes on S, 
normalized so that 

r = t'^ + f^i^^Vw (21) 

where is the stationary Killing field (with unit norm at infinity) and summation over 
/i is understood. (This equation both picks out a particular family of axial Killing fields, 
V?"^), acting in orthogonal planes, and defines the "angular velocity of the horizon", ^^h^ ■) 
Choose Va to be the derivative operator of this solution, so that Va is invariant under 
the isometries generated by Then eq.(|l3D holds, and, in addition, the right side 
now vanishes since £^0 = 0. Let 6(j) be an arbitrary, asymptotically fiat solution of the 
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linearized equations. Then, the fundamental identity eq.(]T^ yields simply 

di5Q) = • 0) 



(22) 



Choose C be an asymptotically flat hypersurface with "interior boundary" S. Integrating 
eq.(0) over C, taking into account eqs. (|16]), (p!?]), and ( pi]) together with the fact that 
C,"" vanishes on S, we obtain 
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f Q = 6S- n^^^6J^,) (23) 



Equation (|23|) corresponds precisely to the first law of black hole mechanics as derived 
by Hamiltonian methods However, eq.(^) has the advantage over this previous 
derivation that the surface term arising from the black hole has now been explicitly 
identified as the variation of the Noether charge of S. 



Equation (23) still is not of the desired form in the sense that the left side of eq.(^) has 
not yet been written as k times the variation of a local, geometrical quantity on S, since Q 
is locally constructed from and its derivatives as well as from the fields appearing in the 
Lagrangian. However, we now will show that the desired form of the first law holds when 
we further restrict attention to the case where 6(f) describes a perturbation to a nearby 
stationary black hole. First, we note that any derivative, Vai...Va„^^ of any Killing field 

can be re-expressed in terms of a linear combination of C,"" and its first derivative, 
Va^fe, with coefficients depending upon the Riemann curvature and its derivatives (see, 
e.g., eq.(C.3.6) of [0). Next, we note that on E we have = and VaC,b = i^^ab, where 
Cab denotes the bi-normal to S. Now define the {n — 2)-form Q on S by the following 
algorithm: Express Q in terms of and Va^b by eliminating the higher derivatives of 
as described above. Then set = and replace Va^b by eab- Since any reference to 
has been eliminated, we see that Q is locally constructed out of the fields appearing in L. 
Furthermore, since Q on S is determined by a well defined algorithm whose only input is 
a Lagrangian L which is invariant under diffeomorphisms of the dynamical fields, 0, (see 
eq.(^ above) it follows that Q is similarly invariant under spacetime diffeomorphisms 
which map S into itself. Thus, Q is a "local, geometrical quantity" on E. Finally, it is 
worth noting that Q is just the Noether charge {n — 2)- form associated with the Killing 
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field ^'^ = K^", i.e., is the horizon Kilhng field normalized so as to have unit surface 
gravity. 

Now identify the unperturbed and perturbed stationary black hole spacetimes in 
such a way that the Killing horizons of the two spacetimes coincide, and the unit surface 
gravity horizon Killing fields coincide in a neighborhood of the horizons. (That this 
always can be done follows from the general formula for Kruskal-type coordinates given 
in Note, however, that for a perturbation which changes the surface gravity, we 

cannot identify the two spacetimes so that the two horizon Killing fields coincide on the 
horizon when normalized via eg. (pT]) ; in addition, since we take = ^V^"^) = near 
infinity, for a perturbation which changes the requirement that = precludes 

us from choosing even to be proportional to the horizon Killing field near infinity in 
the perturbed spacetime.) Then, on S we have, 

SQ = k5Q (24) 

where k is the surface gravity of the unperturbed black hole. Hence, for perturbations to 
nearby stationary black holes, the first law of black hole mechanics (p3D takes the form 

^6S = 58- n^pSJ^,) (25) 

where the "black hole entropy", 5*, is defined by 

S = 2tt [ Q (26) 

JT, 

Thus, we have established the existence and "local, geometrical character" of the notion 



of black hole entropy, S, in a general theory of gravity [|T8 . 

Our "local, geometrical" formula ( P^D for the entropy of a stationary black hole sug- 
gests the following generalization to the non-stationary case: For an arbitrary cross- 
section, S', of the horizon of a non- stationary black hole, construct Q by exactly the 
same mathematical algorithm as used above for the bifurcation surface, S, of a Killing 
horizon. Then 27r times the integral of Q over S' yields a candidate expression for the 
black hole entropy at "time" S'. The viability of this proposed definition is presently 
under investigation [Q. 
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Note that eq.(P3|) has been derived only for the case of perturbations to nearby sta- 
tionary black holes, even though eq.(^) holds in the more general case of non-stationary 
perturbations. However, since 5^° = and, on S, we have ^'^ = 0, it follows that 
(5[V6^°] = on S. From this, it follows that for non-stationary perturbations, we have, 
on S, 

(5[V[a6]] = l^^(^ab + Wab (27) 

where, as before, eab denotes the binormal to S, and Wab is purely "normal-tangential", 
i.e., it vanishes when both of its indices are projected into S or both projected normal 
to S. It then follows from the existence of a reflection isometry about S (see lemma 2.3 



of [|l^]) that the Wab-term makes no contribution to the variation of Q. It then can be 
seen that for sufficiently "low derivative" theories where Q depends only upon and 
its first antisymmetrized derivative (as occurs, in particular, in general relativity and in 
(1 + l)-dimensional theories), eg. (|2^) continues to hold for non-stationary perturbations. 
Thus, in such theories, the first law of black hole mechanics continues to hold in the form 
(|25|) , with 5* given by (|26D . The nature of the first law for non- stationary perturbations 
in more general theories is presently under investigation 

The fact that for stationary black holes, S is just 2tt times the Noether charge of 
the horizon Killing field (normalized to have unit surface gravity) implies that for an 
initially stationary black hole which undergoes a dynamical process and later "settles 
down" to a stationary final state, the net change in black hole entropy is just the total 
fiux through the horizon of Noether current associated with a suitable "time translation" 
on the horizon. This suggests a possible relationship between the validity of the second 
law of black hole mechanics in a theory and positive energy properties of that theory. 
It also suggests some possible approaches toward establishing (or disproving) the second 
law in general theories of gravity. These issues also are under investigation 

Finally, we consider the relationship between the results of this paper and the formula 
for black hole entropy obtained via the "Euclidean approach" in the manner first given 
in 0. We begin by noting that since = 0, the Noether current (|^) associated with 
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the horizon Kilhng field, of a stationary black hole is simply 



j = -e • L (28) 

Let C be an asymptotically fiat hypersurface with "interior boundary" S. Integrating 
eq.(^) over C and taking into account eqs.(|^), (0), (0), and (|20|), we obtain, 

J^Q = -J^^.L- Jj-B (29) 

Now, the "Euclidean action" , /, corresponds to 

-[/ e-L+ / t-B] (30) 

K Jc Joo 



More precisely, in the static case, the right side of eq. (|30D equals what would be obtained 
by integrating the suitably analytically continued Lagrangian, L' = L + dB, over a 
"Euclidean section", constructed by replacing the Killing parameter, t, by r = it, and 



then periodically identifying r with period 27t/k (see, e.g., [19| for further details). In 
the stationary but non-static case, there is no such thing as a "Euclidean section", but 
the right side of eq. ( ^OD corresponds to what researchers mean by the "Euclidean action" 
in that case. Thus, we obtain the following formula for /, 

^I = S-nPj^^^- J^Q (31) 

Now, in the Euclidean approach, is identified as the thermodynamic potential of the 
black hole [0] . This leads immediately to the the following formula for black hole entropy. 



s 



Q (32) 



which agrees with eq. (|26D . Thus, we have shown that the "Euclidean procedure" for 
obtaining black hole entropy gives the same result as obtained by our method. 

This research was supported in part by National Science Foundation grant PHY- 
9220644 to the University of Chicago. 
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